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Abstract 

We discuss the role of the Zweig rule (ZR) violation in the scalar channel for the 
determination of low-energy constants and condensates arising in the effective chiral 
Lagrangian of QCD. The analysis of the Goldstone boson masses and decay con- 
stants shows that the three-flavor condensate and some low-energy constants are 
very sensitive to the value of the ZR violating constant Lq. A similar study is per- 
formed in the case of the decay constants. A chiral sum rule based on experimental 
data in the ++ channel is used to constrain Lq, indicating a significant decrease 
between the two- and the three-flavor condensates. The analysis of the scalar form 
factors of the pion at zero momentum suggests that the pseudoscalar decay constant 
could also be suppressed from Nf = 2 to 3. 



1 Introduction 



The low-energy constants (LEC's) of the effective chiral Lagrangian of QCD 
[1] are quantities of great theoretical interest, since they reflect the way chiral 
symmetry is spontaneously broken. However, their determination remains a 
particularly awkward problem. In most cases [1,3-6], their values have been 
inferred from observables for the pseudoscalar mesons, with the help of two 
assumptions: (1) the quark condensate is the dominant order parameter to 
describe the Spontaneous Breakdown of Chiral Symmetry (SB^S) [1], and (2) 
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the pattern of SB^S agrees correctly with a large- N c description of QCD [7], 
in which quantum fluctuations are treated as small perturbations. 

If we admit both assumptions, the SU(2)xSU(2) quark condensate S(2) = 
— \im mu ^ md ^ {uu} should not depend much on the mass of the strange quark. 
We could then set the latter to zero with no major effect on the quark conden- 
sate: S(2) ~ lim ms ^o ^(2) = S(3). We end up with only one large condensate 
for SU(2)xSU(2) and SU(3)xSU(3) chiral limits, which is not very sensitive 
to qq fluctuations. The LEC's suppressed by the Zweig rule, L 4 and L 6 , are 
consistently supposed to be very small when considered at a typical hadronic 
scale. 

However, several arguments may be raised against this "mean- field approxi- 
mation" of SB^S, in which the Zweig rule applies and the chiral structure of 
QCD vacuum is more or less independent of the number of massless quarks. On 
the one hand, the scalar sector ++ does not comply with large iV c -predictions 
[8], and some lattice simulations with dynamical fermions suggest a strong 
iV/- sensitivity of SB^S signals [9]. On the other hand, the behaviour of the 
perturbative QCD /5-function indicates that chiral symmetry should be re- 
stored for large enough values of Nf. In the vicinity ot the corresponding 
critical point, chiral order parameters should strongly vary with Nf. Various 
approaches, based on the investigation of the QCD conformal window [10], 
gap equations [11], or the instanton liquid model [12], have been proposed 
to investigate the variations of chiral order parameters with the number of 
massless flavors and to determine the critical value of Nf for the restoration 
of chiral symmetry. 

In Ref. [13], the ^/-sensitivity of chiral order parameters has been investigated 
without relying on perturbative methods, but rather by exploiting particular 
properties of vector-like gauge theories. The mechanism of SB^S is indeed 
related to the dynamics of the lowest eigenvalues of the Dirac operator: Jj) = 
luid^ + iGfj), considered on an Euclidean torus [14-16]. Two main chiral order 
parameters can be expressed in this framework. The quark condensate £ is 
related to the average density of eigenvalues around zero [14] and the pion 
decay constant in the chiral limit F can be interpreted as a conductivity [16]. 
The paramagnetic behaviour of Dirac eigenvalues leads to a suppression of 
both order parameters when the number of flavors increases: 

F 2 (iV / + l)<F 2 (iV / ), Z(Nf + 1) < V(N f ) (1) 

This sensitivity of chiral order parameters to light-quark loops is suppressed in 
the large- A^ c limit and is considered as weak for QCD according to the second 
hypothesis of the Standard framework. However, the iV/-dependence of chiral 
order parameters can be measured by correlators that violate the Zweig rule 
in the scalar (vacuum) channel. For instance, the difference E(2) — E(3) (and 
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the LEC Lq) is related to the correlator (uu ss) [17,18] (this correlator can 
be interpreted as fluctuations of the density of Dirac eigenvalues [13]). The 
large ZR violations observed in the ++ channel could therefore support a 
swift evolution in the chiral structure of the vacuum from Nf = 2 to Nf = 3. 
The quantum fluctuations of qq pairs would then play an essential role in the 
low-energy dynamics of QCD. 

Hence, it is worth reconsidering the determination of LEC's without suppos- 
ing (1) the dominance of the quark condensate and (2) the suppression of 
quantum fluctuations. This determination starts with quark mass expansion 
of measured observables such as F 2 M 2 or F\M\, using Chiral Perturbation 
Theory ( X PT) [1]: 

F 2 P M 2 P = m q Y> P + m\\a P + 6 P log(M P )] + F 2 S P , (2) 

where m q denotes formally light quark masses (u, d, s) and the remainder 
Fp5 P is of order rr? q . The coefficients E P , a P , b P are combinations of LEC's. 
The chiral logarithms log(Mp) stem from meson loops. The coefficient of each 
power of m q does not depend on the renormalization scale of the effective 
theory [FpMp and m q are independent of this scale). 

Series like Eq. (2) are assumed to converge on the basis of a genuine dimen- 
sional estimate [19]. The LEC's involved in the coefficients are related to Green 
functions of axial and vector currents, and scalar and pseudoscalar densities. 
The dimensional estimate consists in saturating the correlators by the ex- 
change of resonances with masses of order Aq C d [4]. We obtain coefficients of 
order ~ 1/Aq CD for the power m™. The quark mass expansion would therefore 
lead to (convergent) series in powers of m g /AqcD <^ 1- 

Notice that this genuine estimate cannot be applied to the linear term Ep 
corresponding to the quark condensate (there is no colored physical state to 
saturate (qq))- Moreover, the convergence of the whole series does not imply 
that the linear term Ep is dominant with respect to the quadratic term. In 
this article, we will precisely address (1) the possibility of such a competition 
between the first two orders in the quark mass expansions, and (2) the implica- 
tions of large values for the ZR-suppressed constants L 4 and L 6 , in particular 
for the determination of LEC's. 

Unfortunately, the masses and decay constants of the Goldstone bosons do not 
provide enough information to estimate the actual size of quantum fluctuations 
in QCD. To reach this goal, Refs. [17,18] have proposed a sum rule to estimate 
Lq [or E(2) — E(3)] from experimental data in the scalar channel. Starting with 
Standard assumptions (two- and three-flavor condensates of large and similar 
sizes), Ref. [17] ended up with a ratio E(3)/E(2) ~ 1/2 at the Standard 0(p 4 ) 
order, whereas Ref. [18] confirmed a large decrease of the quark condensate 
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when Standard 0(p 6 ) contributions were taken into account. Even though 
these results suggest a significant variation in the pattern of SB^S from Nf = 2 
to Nf = 3, it is seems necessary to reevaluate this sum rule without any 
supposition about the size of the condensates. This analysis will be performed 
in the second part of this article. 

We will follow mainly the line of Ref. [20], which can be considered as an 
orientation guide to this article. The first part is devoted to the determination 
of the LEC's from the pseudoscalar spectrum. Sec. 2 considers the role played 
by Lq for the Goldstone boson masses and the quark condensates, whereas 
the decay constants and L4 are treated in Sec. 3. Sec. 4 deals mainly with 
the implication of ZR violation in the ++ channel for the determination of 
LEC's. The second part of this article focuses on the estimate of L 6 from data 
in the scalar sector. Sec. 5 introduces the sum rule for L 6 , and sketches the 
Operator Product Expansion of the involved Green function. In Sec. 6, we 
estimate this sum rule, with a special emphasis on the the scalar form factors 
of the pion and the kaon. In Sec. 7, we present the results obtained for the 
quark condensates and LEC's from the sum rule, and we discuss two other 
quantities related to the pion scalar form factors: the slope of the strange form 
factor and the scalar radius of the pion. Sec. 8 sums up the main results of 
the article. App. A collects the expansions of pseudoscalar masses and decay 
constants in powers of quark masses. App. B deals with the Operator Product 
Expansion of the correlator (uu ss). App. C provides logarithmic derivatives 
of the pseudoscalar masses with respect to the quark masses. 



2 Constraints from the pseudoscalar meson masses 

2.1 Role of L 6 

Let us first study the pseudoscalar masses M n , Mk, M v , starting from their 
expansion at the Standard 0{p A ) order, Eqs. (10.7) in Ref. [1]. We reexpress 
them as: 

F 2 M 2 = 2mS(3) + 2m(m s + 2m)Z s + Am 2 A + 4m 2 B 2 L + F 2 5 w , (3) 
F 2 K M 2 K = (m s + m) E (3) + (m s + m) (m s + 2m) Z s (4) 
+ (m s + m) 2 A + m(m s + 7ti)BqL + F^5 K , 

where m = {m u + md)/2 and Z s and A are scale-independent constants, 
containing respectively the LEC's L 6 (/x) and L 8 (fj,), 
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with i?o = S(3)/F 2 and Fo = -F(3). The remaining 0(p 4 ) chiral logarithms 
are contained in LF£] 
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There is a similar equation for F 2 M 2 : 



F 2 M 2 = -(2m, + m)S + -(2m a + m)(m, + 2m) 



(8) 



+^(2m 2 + m 2 )A + ^(m s 



with the scale-independent constant Z p = 16BqL 7 . A factor B is included in 
the expression of A, Z$ and Zp in terms of Lj =6 7 8 , so that they do not diverge 
in the limit S(3) — > 0. The corresponding equations for the pseudoscalar decay 
constants Fp (P— it, K, 77) will be treated in Sec. 3.1. 

We take FpMp and Fp as independent observables, in order to separate in 
a straightforward way the "mass" constants Lq, L7, Lg from L4, L5 that ap- 
pear only in the expansion of decay constants Fp. There is a second argument 
supporting the choice of Fp and FpMp as independent observables of the 
pseudoscalar spectrum. We expand observables in powers of quark masses, 
supposing a good convergence of the series. We have sketched in the intro- 
duction how a naive dimensional estimate justifies this assumption : LEC's 
are related to QCD correlation functions, which can be saturated by massive 
resonances, leading to series in (m q / Aqcd)- We should therefore expect good 
convergence properties for "primary" observables obtained directly from the 
low-energy behavior of QCD correlation functions, like Fp and FpMp. For 
such quantities, the higher-order remainders should thus remain small. On 
the other hand, we have to be careful when we deal with "secondary" quanti- 
ties combining "primary" observables. The higher-order remainders may then 
have a larger influence. In particular, ratios of "primary" observables (like 
Mp = FpMp/ Fp) might be dangerous if higher-order terms turned out to be 
sizable [leading to untrustworthy approximations like 1/(1 + x) ~ 1 — x with 
a large x\. 



2 In this article, we use the following values of masses and decay constants: M n 
135 MeV, M K = 495 MeV, M v = 547 MeV, F n = 92 A MeV, F K /F n = 1.22. 
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In Eqs. (3), (4) and (8), all terms linear and quadratic in quark masses are 
shown. The remaining contributions, of order 0(m;L ark ) and higher, are col- 
lected in the remainders Sp. We can consider that the latter are given to us, 
so that Eqs. (3), (4) and (8) can be seen as algebraic identities relating the 
3-flavor condensate S(3) the quark mass ratio r = m s /m, and the LEC's 
Fq, Lq(h), Lj and Lg(yu). The three-flavor quark condensate is measured in 
physical units, using the Gell-Mann-Oakes ratio: X(3) = 2mS(3)/(F 7r M 7r ) 2 
[21]. 

We are going to assume that the remainders 8 P are small (S P <C Mp), and 
investigate then the consequences of Eqs. (3) and (4) for the values of LEC's. 
Before starting, we should comment the status of Eqs. (3), (4) and (8) with 
respect to Chiral Perturbation Theory (%PT). Even if we imposed Sp = 0, 
we would not work in the frame of one-loop Standard %PT [1]: we do not 
suppose that the condensate S(3) is dominant in these equations, we do not 
treat 1 — X(3) as a small expansion parameter, and accordingly, we do not 
replace (for instance) 2itlBq by M 2 in higher-order terms. However, we are not 
following Generalized either [22] , since Bq is not treated as an expansion 
parameter: even with Sp = 0, Eqs. (3), (4) and (8) exceed the Generalized tree 
level, since these equations include chiral logarithms. 

It is useful to rewrite Eqs. (3) and (4) as: 



2^ x / x <? n Am 2 B 2 rL 

[£(3) + (2m + m s )Z s ] = 1 - e(r) - F2 ° - - S, (9) 



F 2 M 2L w ■ v FlMl r-1 



Am 2 A . Am 2 B 2 L . 

T^TTTZ , + ( 10 ) 



with 



F 2 M 2 v ' F 2 M 2 r - 1 

^ TV 1V± TT ^ TV 1V± TT 



J 2 -r _ n (F K M K ^ 2 



S and 8' are linear combinations of the remainders 8 n and Sk'- 



* = ^l^-| l + r^r)^, (12) 

8-K 
TV 

For large r, we expect 8' <C 8 ~ 8^/M 2 . Similarly to Ref. [6], we consider as 
parameters F Q = \im mrns _> F n [i.e. L 4 (//)], the ZR violating constant Lq(h) 
and the quark mass ratio r = m s /m. Eq. (9) ends up with a non-perturbative 
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formula (no expansion) for the three-flavor Gell-Mann-Oakes-Renner ratio 
X(3): 



where k contains L 6 (/i): 



'F W M W \ 2 



« = 64(r + 2)1-^) {LM (15) 



o 



' M K 2 M„\ rL 
log h - log — - + 



256tt 2 V & fi 9 & fj, J 16(r-l)(r + 2) 



Eq. (14) is an exact identity, useful if the remainder 8 in Eq. (12) is small, i.e. 
if the expansion of QCD correlators in powers of the quark masses m„, ma, 
m s is globally convergent. It means that 5p <C Mp in Eqs. (3) and (4), but 
the linear term in these equations (related to the condensate) does not need 
to dominate. 

k describes quantum fluctuations of the condensate, and actually k = 0(1/N C ). 
L 6 (/jl) has to be fixed carefully to keep k small, k is equal to zero for 10 3 • L 6 = 
—0.26 at the scale \x = M p , which is close to the value usually claimed in 
Standard analysis [3]. In this case, Eq. (14) yields X(3) near 1, unless 
the quark mass ratio r decreases significantly, leading to e — > 1. This effect is 
well-known in G%PT [22]: the minimal value of r (corresponding to e = 1 and 
X(3) = 0) is fx = 2(F K M K )/(F 7T M 7T ) - 1 ~ 8. Notice that for these very small 
values of r, the combination of higher-order remainders 8 cannot be neglected 
any more in Eq. (14). 

But quantum fluctuations can modify this picture: the number before the 
curly brackets in Eq. (15) is very large (~ 5340 for r = 26 and F = 85 MeV). 
Hence, even a small positive value of L 6 (M p ) can lead to a strong suppression 
of X(3), whatever the value of r = m s /m. This effect can be seen on Fig. 1, 
where X(3) is plotted as a function of L 6 (M P ) for r = 20, r = 25, r = 30 and 
F = 85 and 75 MeV. The decrease of X(3) is slightly steeper for lower values 
ofF . 

Once X(3) is known, Eq. (10) leads to L 8 (pi): 

r / x L F* e + 8 1 1 f, Ml 2, Af 2 l , x 
LM = — + F^Ml [X(3W + 512^ l l0g # + 3 lQg #} • (16) 
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Fig. 1. X(2) (crosses) and X(3) (no symbol) as functions of Lq(M p ) ■ 10 3 and 
r = m s /m (solid line: r = 20, dashed line: r = 25, dotted line: r = 30) for ^0=85 
MeV (upper plot) and 75 MeV (lower plot). 

This constant depends on L 6 only through X(3). Notice that this dependence 
is smaller when F decreases (L 8 depends actually on L 6 through B ). 

The ZR violating constant L 7 can be obtained from F^M^ Eq. (8): 
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where F v will be discussed in Sec. 3.1. The pseudoscalar spectrum satisfies with 
a good accuracy the relation 3F 2 M 2 + F 2 M 2 = 4F^M^, which reduces at the 
leading order to the Gell-Mann-Okubo formula [23]. This relation leads to a 
strong correlation between A and Zp\ A + 2Zp ~ 0. This correlation can also 
be seen in Standard %PT between Lj and L$, and remains to be explained in 
both frameworks. No obvious reason forces this particular combination of two 
low-energy constants to be much smaller than the typical size of the effective 
constants. 



2.2 Paramagnetic inequality for E 



In Ref. [13], qq fluctuations were shown to increase the two-flavor condensate 
S(2) = — \im m ^o(uu) \ ms physical with respect to 11(3), so that X(2) > X(3). 
The two-flavor quark condensate can be obtained through the limit: 

S (2) = lim (F ; M?r)2 = E(3) + m s Z s \ m=0 + S 2 , (18) 



keeping m s fixed. We have the quantities: 82 
Z s + B$AZ S , with: 



lim m ^ F%8ir/ (2m) and Z' 
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, m k 2 , Ml 
log + - ■ log -=2- 



(19) 



and Mp = lim m ^ Mp. The effect of AZ S is very smallQ. AZ S should be 
compared to the logarithmic terms included in Z s , Eq. (5), at a typical scale 
/i ~ M p . AZ S reaches hardly 10% of this logarithmic piece. 

Once Z s is eliminated from Eqs. (9) and (18), we obtain the two-flavor Gell- 
Mann-Oakes-Renner ratio X(2) = 2mE(2)/(F 7r M^) 2 : 



X(2) = [l- e 1— + — JTO) 
{F n M,Y 



(20) 



2F 4 



-x(sy 
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(r - l)(r + 2) 



L - rAZ s 



+ 8 



with: 



It can be evaluated following the procedure of Sec. A. 3. 
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In the expression of 5x, the remainders 5p/Mp are suppressed by a factor 
m/m s : this suppression is obvious for 5^ [e = 0(l/r)], whereas the operator 
applied to ^ cancels the terms of order 0(mm 2 s ). For r > 20, we expect thus 
\5x\ ~ The dependence of X(2) on L 6 is completely hidden in X(3), and 
therefore marginal, as shown in Fig. 1. 

The paramagnetic inequality X(2) > X(3) constrains the maximal value 
reached by X(3) = X(3)\ max . If we neglect <5£(2) in Eq. (18), the inequal- 
ity can be translated into a lower bound for L 6 : 



5127T 2 



Ml 
Ai 2 
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(23) 



Fig. 1 shows clearly the lower bound: Lq(M p ) > —0.21 • 10 3 . 

X{2) is loosely related to X{3), but it is very strongly correlated with r, 
specially for small values of r. Eq. (20) yields the estimate X(2) ~ [1 — e] • 
r/(r + 2), up to small correcting terms due to X(3). We are going to study 
the effect of these correcting terms. 

If we neglect AZ S , X(2) is a quadratic function of X(3), which is not monotonous 
when X(3) varies from to X(3)| max : it first increases, and then decreases (see 
Fig. 1). The decrease of X{2) for X{3) close to its upper bound is caused by 
the negative term, quadratic in X(3), in Eq. (20). This decrease of X(2) is 
more significant for small F Q , because the factor in front of [X(3)] 2 in Eq. (20) 
becomes larger. 

Therefore, X(2) does not reach its maximum for the paramagnetic bound 
X{2) = X(3)\ max , whereas its minimum is the smallest of the two values 
obtained for X(3) = and X(3) = X(3)| max Q. The dependence on r of the 
minimal value of ^(2) can be guessed rather easily. For large r, the term 
linear in X(3) in Eq. (20) can be neglected: the minimum of X(2) occurs 
for X(3) = X(3)| max . For small r, X{2) and X{3) tend to 0, and the term 
quadratic in X(3) should be small with respect to the linear term. Therefore, 
X(2) reaches its minimum for X(3) = when r is small. 

4 This updates Ref. [13], where the minimum and maximum of X(2) were claimed 
to be obtained for X(3) =0 and X(3) = X(3)| max . 
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Fig. 2. X(2) as a function of r for several values of Fq 

The numerical analysis of Eq. (20), including AZ S , supports this intuitive 
description. In Fig. 2, the variation ranges of X(2) are plotted for several 
values of F . The curve for the minimum of X(2) exhibits a cusp when the 
minimum of X(2) corresponds no more to X(3) = 0, but to X(3)| max . X(2) 
appears to be strongly correlated to r, even though a large value of X{2) ~ 0.9 
can be associated to a broad range of r. 



3 Constraints from the pseudoscalar decay constants 



3.1 Role of L A 



The decomposition used for the Goldstone boson masses can be adapted to 
the decay constants: 



Fl = F 2 + 2m£ + 2(2m + m s )| 
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with the scale- independent constants related to L 4 and L 5 : 
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Eqs. (24) and (25) contain all the terms constant or linear in quark masses in 
the expansion of F 2 and F|-, whereas ep denote remainders of order 0(m 2 uark ). 
There is also a formula for F v , which can be written as: 
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The two scale-independent constants can be extracted from Eqs. (24) and 
(25): 
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(30) 



with: 
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(31) 



where the latter estimate is obtained for F K /F W = 1.22. 



£ (i.e. L 5 ) turns out to depend essentially on X(3) and r, whereas £ (i.e. L 4 ) 
is related to the difference between F(3) and F n . Eq. (30) leads to a quadratic 
equation for [F(3)/F,-] 2 , involving L 4 : 



(f>y- (1 -^>(f>)v,<3Ho, 
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with: 
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(34) 



Eq. (32) has the solution: 



'F(3)\ 2 _ 1 - fj - e + ^(1 - f] - ef - 4AX(3) 



(35) 



Notice that this formula is very close to Eq. (14), that relates X{3) to Lq(h) 
through the parameter k. For r = 25, the factor in front of the curly brackets 
in the definition of A is of order 460, and A vanishes for L 4 (M P ) = —0.51 • 10 -3 . 

Eq. (32) admits a solution only if XX (3) < (1 - fj - ef/A. From Eq. (35), we 
get then a range for F(3): 



1 — f] — e 



< 1 — fj — e. 



(36) 



The variations of F(3) with L±(M p ) are plotted for three values of r and 
X{3) = 0.9 and 0.5 in Fig. 3. When X(3) decreases, the allowed range 
for L 4 broadens. This is due to the definition of L 4 , which relates B L 4 = 
S(3)L 4 /F 2 (3) to the low-energy behavior of a QCD Green function. [F(3)/F 1T ] 2 
starts at 0.9 (for the lowest possible value of L 4 (M p )) and decreases until 0.5. 
F(3) can thus vary from 87 to 65 MeV, depending on the value of L 4 (M P ). 

3.2 Paramagnetic inequality for F 2 

We obtain F(2) by taking the limit: 

F 2 (2) = lim F 2 = F(3) 2 + 2mi| m=0 + e 2 , (37) 
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Fig. 3. F{2) (crosses) and F(3) (no symbol) as functions of L^Mp) ■ 10 3 and 
r = m s /m (solid lines: r = 20, dashed lines: r = 25, dotted lines: r = 30) for 
,X"(3)=0.9 (upper plot) and 0.5 (lower plot). 



keeping m s fixed. We have £2 = linim^o^Tr and £| m =o = £ + B A^, with 



1 1 M| 

32tt 2 g M 2 K 



(38) 
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A£ has a tiny effect on the results, similarly to AZ S for X(3). We get the 
equation: 



F{2) 



1 2 



T I 2 
r + 2| [1 -^ ] + r 



F(3) 



n 2 



X(3)0 



F(3) 



(39) 



with: 



0: 



r + lejr r + 2 e 2 



F 2 

^ 7T 



IF,? 

TV 

"r + 1 



r 



32tt 2 F 2 



Ar 



TV 

r + 2 



r-lj F 2 



K 



lim e 7 

m— >0 



7/ + 



r — 1 



1 M£ 2r 



Ml 

TV 



r — 1 



1 °g^2 L + ( r + 2 ) 1 °g 



Ml 



Ml 



(40) 
(41) 

(42) 



It is interesting to compare the expression of F(2) as a function of F(3) with 
Eq. (20) that relates X{2) to X{3). Even though the equations look rather 
similar, we can notice that Eq. (20) is a quadratic function of X(3) whereas 
Eq. (39) involves [F(3)/F w ] 2 and its inverse [F w /F(3)} 2 . Moreover, Eq. (39) is 
an increasing function of F(3), whereas X(2) is not monotonous with X(3). 
On the other hand, e F suppresses the remainders £p/Fp by a factor m/m s , in 
a comparable way to 5x- 

The paramagnetic inequality F(2) > F{3) is translated into an upper bound 
for F(3): 



F(3) \ 2 < 1 - 77 - e F + ^(1 - 77 - e F f - 40X(3) 



(43) 



A quick estimate shows that the condition X{3) < (1 — 77 — e F ) 2 /(Acf)) is 
satisfied for any r between fi and f 2 . This paramagnetic bound corresponds 
to a lower bound for L 4 (//): 

1 M 2 



The term log(M| : ) is very weakly dependent on r, X(3) and F(3). If we scan 
the acceptable range of variation for these three parameters, we obrain the 
lower bound L 4 (M p ) > —0.37 • 10~ 3 . We notice that the curves for F{2) and 
F(3) cross each other at this value in Fig. 3. 
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F(3) [MeV] 


F{2) [MeV] 


No ZR violation 


r 


min 


max 


min 


max 


-T [Me V J 


A 


1 n 
1U 


bi.fay 


O T O A 

87.24 


81.31 


O T O A 

87.24 


85.57 


0.403 


lb 


no 
63.02 


on 1 o 
89.12 


CO OA 

82.34 


on 1 O 

89.12 


86.1 / 


0. /ol 


on 


63.63 


89.99 


83.02 


89.99 


86.71 


0.860 


25 


63.99 


90.50 


83.38 


90.50 


87.08 


0.905 


30 


64.23 


90.83 


83.43 


90.83 


87.35 


0.928 


35 


64.39 


91.06 


83.14 


91.06 


87.55 


0.940 



Table 1 

On the left hand-side, bounds for F(3). In the middle, corresponding bounds for 
F{2). On the right hand-side, values of F(3) and X(3) saturating both paramagnetic 
inequalities: F{2) = F(3) and X(2) = X(3). 



Since > 0, Eqs. (36) and (43) lead to the range: 
1-fj-e ^ F(3) 2 



< 



F 2 

^ TV 



< 1 — rj — max(£p, e). 



(45) 



The bounds on F(3) are indicated in Tab. 1 (neglecting the remainders e 
and £p). Table 1 collects for several values of r the corresponding bounds for 
F(2), obtained using Eq. (39). It gives also the values of X(3) and F(3) that 
saturate both paramagnetic bounds: F(2) = F(3) and X(2) = X(3). In this 
case, the Zweig rule would be violated neither for the masses nor for the decay- 
constants. 



4 Sensitivity of low-energy constants to ZR violation 



The equations (16), (17), (29) and (30) can be used to extract the LEC's 
^=4,5,7,8 as functions of r, F et X(3), or equivalently, of L 6 , F = F(3) and 
r using Eq. (14). Results are shown in Tab. 2 as a function of L 6 , for 2 values 
of F and 3 values of r. 

This table has been obtained by neglecting the higher-order terms 5p and ep, 
which start at the next-to-next-to- leading order (NNLO). If the size of these 
remainders is large, the values collected in these tables should be clearly mod- 
ified. If we keep considering the low-energy constants as functions of L 6 (M p ) 
and we change the relative size of the NNLO remainders within a range of 5%, 
the corresponding variations of X(3) are of order ±0.02. The impact on the 
other LEC's is larger for smaller X(3) (of order ±0.3 • 10~ 3 for X(3) ~ 0.8, 
r ~ 25). 
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r = 


20 






r = 


25 






r = 


30 







4 


5 


7 


8 


4 


5 


7 


8 


4 


5 


7 


8 


-0.2 


-0.284 


2.410 


-1.259 


2.624 


-0.282 


1.603 


-0.503 


0.994 


-0.284 


1.130 


-0.185 


0.298 


-0.1 


-0.264 


2.628 


-1.452 


3.044 


-0.261 


1.812 


-0.604 


1.224 


-0.262 


1.332 


-0.233 


0.416 


0. 


-0.246 


2.824 


-1.636 


3.445 


-0.242 


1.993 


-0.699 


1.440 


-0.243 


1.501 


-0.276 


0.525 


0.2 


-0.215 


3.167 


-1.986 


4.207 


-0.210 


2.306 


-0.880 


1.848 


-0.211 


1.790 


-0.360 


0.730 


0.4 


-0.188 


3.466 


-2.318 


4.930 


-0.184 


2.572 


-1.050 


2.232 


-0.184 


2.036 


-0.439 


0.925 


1. 


-0.121 


4.209 


-3.256 


6.966 


-0.116 


3.228 


-1.532 


3.314 


-0.117 


2.631 


-0.664 


1.471 


-0.2 


0.068 


1.833 


-0.816 


1.653 


-0.013 


1.174 


-0.325 


0.586 


-0.065 


0.782 


-0.117 


0.124 


-0.1 


0.133 


2.088 


-0.999 


2.056 


0.053 


1.414 


-0.420 


0.804 


0.002 


1.014 


-0.161 


0.236 


0. 


0.190 


2.310 


-1.176 


2.442 


0.108 


1.616 


-0.509 


1.008 


0.055 


1.203 


-0.202 


0.340 


0.2 


0.287 


2.684 


-1.503 


3.156 


0.200 


1.953 


-0.678 


1.391 


0.144 


1.513 


-0.290 


0.533 


0.4 


0.368 


3.002 


-1.814 


3.833 


0.277 


2.234 


-0.837 


1.751 


0.217 


1.770 


-0.354 


0.715 


1. 


0.567 


3.781 


-2.696 


5.751 


0.463 


2.915 


-1.291 


2.774 


0.392 


2.384 


-0.565 


1.232 



Table 2. Low- 
part). 



■energy constants Lj(M p ) • 10 3 as functions of Lq{M p ) ■ 10 3 and r = m s /m, for Fq=85 MeV (upper part) and 75 MeV (lower 



The authors of Ref. [6] have estimated the NNLO remainders in the Standard 
Framework. The authors assume first L 4 (M p ) = L 6 (M p ) = and r = 24, 
estimate 0(p 6 ) counterterms (Standard counting) through a saturation of the 
associated correlators by resonances, and fit globally on the available data 
(masses, decay constants, form factors). For the decay constants, the obtained 
NNLO remainders ep are less than 5%. The situation is less clear for the 
masses, due to a bad convergence of the series. For instance, Ref. [6] has 
obtained the decomposition: M 2 /(M 2 ) phys = 0.746 + 0.007 + 0.247, where 
the three terms correspond respectively to the leading 0(m), next-to-leading 
0(m 2 ) and next-to-next-to-leading 0(m 3 ) orders. Ref. [6] suggested that a 
variation of L 4 (M P ) and Lq(M p ) could make the NNLO contribution decrease, 
but a competition occurs then between the 0{m 2 ) term and the leading-order 
term. 

Such a competition could be understood from our analysis as a consequence of 
the suppression of the three-flavor condensate X(3). It would be of great inter- 
est to proceed to the same analysis as in Ref. [6], and to allow a competition 
between the terms linear and quadratic in quark masses. This might improve 
the convergence of the expansion of Goldstone boson observables in powers 
of quark masses. Even if the three-flavor condensate X(3) is suppressed (first 
term of the series for pseudoscalar masses) , we expect a global convergence of 
the series, i.e. small higher-order remainders 5p and ep. 

Standard values of the LEC's at order 0(p 4 ) can be found in Ref. [3], and 
were derived with the supposition that the ZR violating LEC's L Q {ji) and 
L^n) were suppressed at the scale fi = M v . The following values have been 
obtained: L 4 (M P ) ■ 10 3 = -0.3 ± 0.5, L 5 (M P ) ■ 10 3 = 1.4 ± 0.5, L 6 (M P ) ■ 10 3 = 
-0.2±0.3, L 7 (M P )-10 3 = -0.4±0.2, L 8 (M P )-10 3 = 0.9±0.3. These values are 
compatible with our analysis: it can be seen on the first lines of Table 2, with 
L 6 (M p ) ~ -0.2 • 10~ 3 , X{2) ~ X{3) ~ 0.9, r = 25, F = 85 MeV. The values 
of L 4 (M p ) and L 6 (M p ) correspond also to the lower bounds derived from the 
saturation of the paramagnetic inequalities for F 2 and X: X{2) = X(3) and 
F(2) = F(3). 

We see that L 4 is weakly sensitive to L 6 , in agreement with Eqs. (27) and 
(30). For large r, and F close to F n , 1 — fj(r) — F^/F 2 is nearly vanishing, 
so that L 4 is mainly fixed by the last term in Eq. (30) with no dependence 
on X(3). On the contrary, L5 is clearly dependent on the value of Lq. We can 
look at Eq. (29) to understand this phenomenon : fj(r) may be small, but it 
never vanishes. The first term in Eq. (29) leads therefore to large values of L 5 
when X(3) decreases, whatever values we choose for r and F (this is related 
to the definition of L 5 ). The Gell-Mann-Okubo formula Eq. (17) correlates 
strongly L 7 and L 8 , leading to L 7 ~ —L 8 /2. 

From this analysis of the pseudoscalar masses and decay constants, we cannot 
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conclude whether qq fluctuations have important effects on the chiral structure 
of QCD vacuum. However, the values of LEC's are extremely sensitive to the 
value of Lq(M p ). A small shift of L 6 towards positive values would immediately 
split X(2) and X(3), and increase strongly L 8 (M p ) and L 5 (M p ). We will now 
use additional information from experimental data in the scalar sector, in 
order to constrain L 6 (M p ). 



5 Sum rule for X(2) - X(3) 



5. 1 Correlator of two scalar densities 



We introduce the correlator [17,18]: 

U{p2) = i WW^bl * x ^ <°l T ^) ^(0)}|0>, (46) 

7T K 



that is invariant under the QCD renormalization group, and violates the Zweig 
rule in the ++ channel. For m s ^ 0, II is a SUl(2) ® SU#(2) order parame- 
ter, related to the derivative of S(2) with respect to m s : mm s <9X(2) / 'dm s = 

We can use the relation Eq. (20) between r, X(3) and X(2) to compute 
<9£(2)/<9m s . Eq. (5) can be used to compute (dZ s /dm s ) m=0 . This leads to 
an equation involving Z s and 11(0): 



X(2)-X(3) = 2 -^Z S + ^-hm^ (47) 



m=0 F 2 M 2 



m 



2ML ff) , r[X(3)] 2 F 2 M 2 / 2 



F 2 w 32tt 2 F n 4 V 9 




F 2 S^ 

mri v '"(^jSoIt' (48) 



IT I 



with the logarithmic derivatives: Ap = m s • 9(log Mp)/dm s . This equation 
contains the NNLO remainder -F^^ in the quark mass expansion of F 2 M 2 . 
Its leading term is 0(mm 2 ), so that the last term in Eq. (48) should be of 
order ~ (-5 n /M 2 ). 

Lq (or the difference X(2) — X(3)) measures the violation of the Zweig rule 
in the scalar channel. We are going to exploit experimental information about 
this violation and to evaluate 11(0) through the sum rule: 
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Fig. 4. Contour of the integral in the sum rule for the correlator (uu ss). 
The three terms will be estimated in different ways: 

• For < y/s < y/si ~ 1.2 GeV, the spectral function Im II is obtained 
by solving Omnes-Muskhelishvili equations for two coupled channels, using 
several T-matrix models in the scalar sector. 

• For y/s~[ < y/s < y/so ~ 1.5 GeV, the spectral function under s\ is exploited 
through another sum rule in order to bound the contribution of the integral. 

• For \s\ = s , we estimate the integral through Operator Product Expansion 
(OPE). 



5.2 Asymptotic behavior 



IT can be expanded using OPE: 



n(p2) = ij^^J^^" W T lM x )Ss(0)]\0) (50) 

p-»i^S(74^ c< " ,(i)(0|O " |0> ' (51) 

with P 2 = —p 2 , n a renormalization scale, t = fi 2 /P 2 , and O n a combination 
of n-dimensional operators. II transforms chirally as {uu){ss) and we take the 
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chiral limit m — > 0. Hence, the lowest-dimension operator is O4 = m s uu, and 
the contributing diagrams include at least two gluonic lines [17]. 




Fig. 5. Feynman diagrams contributing to m s {uu) in OPE of II (lowest order in a s ). 
The white circle is the scalar source uu, the black one ss. 

We will work in dimensional regularization (d — 4 — 2u). In the class of 
t'Hooft's gauges, the gluon propagator reads: 

(52) 

with £ a free real parameter. The Wilson coefficient of m s uu (at the leading 
order) is obtained by adding 6 two-loop integrals. It is easy to see that the 
contributions of £ and £ 2 cancel in this sum of integrals. Hence, we check 
that the Wilson coefficient of m s uu at the leading order is independent of the 
chosen gauge, in the classe of t'Hooft's gauges. 

We want the large-P 2 limit of integrals like: 
1 f d 4 q d 4 k 

g 4 s fi 4uJ m s (uu} — J -^-^-^-^V(p 2 ,q 2 ,k 2 ,p- q,p- k,q- k,m 2 s ) (53) 
1 

X 

[(p - q) 2 - ml] n *[q 2 - m 2 ,] 713 

1 

[(k + p) 2 - m 2 ]«4 [(fc + g)2 _ m 2]n 5 [p _ m 2]n 6 ' 

where V is a polynomial of degree 2. m corresponds at the same time to 
m = m u = rrid for fermion propagators in the loop of u — d quarks, and to a 
fictitious mass to regularize infrared gluonic divergences (we take at the end 
the limit m — > 0). 
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Using identities like 2(k-q) = [(k + q) 2 + m 2 ] — [k 2 + m 2 ] —q 2 , we can reexpress 
the sum in terms of: 



1 r/r i r i \ 1 f d 4 q d 4 k . . 



p 2z/ ° p 2j/ ° 7 [g 2 - m 2 ,]^ [k 2 - m 2 ]^ 

1 



x 



[(k + q) 2 - ml]^[(p - q) 2 - rag]" 4 [(A; + p) 2 -m 2 ]^' 



with mi = m 4 = m and m 2 = = m 5 = m s . These integrals are formally 
identical to the ones arising for two-loop self-energies. The behavior of such 
integrals at large external momentum has already been studied. The basic 
idea consists in following the flux of the large external momentum through 
the Feynman diagrams, in order to Taylor expand the propagators [24]. This 
procedure, based on the asymptotic expansion theorem [25], is sketched in 
App. B. 

Rather lengthy computations lead to the first term arising in the OPE of the 
correlator. Some integrals contain poles in 1/cu, but these divergences cancel 
when all the contributions are summed (this cancellation is a non-trivial check 
of the procedure). The first term in OPE is: 



,— — 2- lim / d A x e %p ' x (0\T{muu(x) m s ss(0)}|0) (55) 

7T K 



P^oo P 2 V 7T / MIM 2 R 



18[1-2C(3)] (a s \ 2 m 2 



m(uu). 



The involved condensate should be the two-flavor one (m — >• 0, m s ^ 0). 

5.3 Contribution for s < si : pion and kaon scalar form factors 

5.3.1 Omnes-Muskhelishvili equations 
In order to compute the integral: 



si 

J = — fdslm U(s) - ( 1 - — V (56) 

7T J S \ So J 



we have to know Im II between and s± (\fs 1 ~ 1.2 Gev). The procedure is 
explained in detail in Refs. [17,18], and we shall merely sketch its main features 
for completeness. In the range of energy between and s±, the im- and KK- 
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channels should dominate the spectral function [17,26,27]. If we denote these 
channels respectively 1 and 2, the spectral function is: 

Im U(s) = E ]l^^[n l F i (s)}[n l G*( S )}e(s - AM?), (57) 



with the scalar form factors for the pion and the kaon: 



F(s) 



(0\uu\tt7t) 
K (0\uu\KK) 



( 



G(s) 



(0\ss\tyty) 



{(0\ss\KK)) 



(58) 



with Mi = M w and M 2 = M K . n x = y 3/2 and n 2 = V% are numerical factors 
related to the normalization of the states and \KK). 

The form factors are analytic functions in the complex plane, with the ex- 
ception of a right cut along the real axis. They should have the asymptotic 
behaviour Fi(s) ~ 1/s when s — > oo, and verify a dispersion relation with 
no subtraction. Obviously, when s increases, new channels open, and the two- 
channel approximation is no more sufficient. But we need F and G for s < s\, 
and we are not interested in the behaviour of the spectral function at much 
higher energies. We can therefore suppose that the two-channel approximation 
is valid for any energies, with the discontinuity along the cut: 



Sij = Sij + 2ia] l2 T ij( y} ,2 d(s - AMf)d(s - AM]) , (59) 

n fl _ am 2 
Im Fi(s) = ^2T*j(s)<jj(s)Fj(s)6(s — AM 2 ), a i = J* — , (60) 

3=1 V S 

and we will suppose in addition that the two-channel T-matrix model impose 
the correct asymptotic behaviour for F and G. 

Under these assumptions, F and G satisfy separately a set of coupled Omnes- 
Muskhelishvili equations [17,26,28,29]: 

1 n °r 1 / s' - AM 2 

m = -ll J d^T; j (s')^-- r ^9(s' -AM^s'), (61) 



with the condition that the T matrix leads to the expected decrease of the 
form factors for s — > oo. Ref. [17] has proved a condition of existence and 
unicity for the solution of Eq. (61): A(oo) — A(AM 2 ) = 2tt, where A(s) is the 
sum of the nn and KK phase shifts. In that case, the set of linear equations 
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admits a unique solution, once the values at a given energy are fixed [18]. All 
the solutions are thus linear combinations of a basis, for instance the solutions 
A(s) and B(s) such as: 1(0) = (J) and B(0) = (J). F and G can therefore be 
written as: 

F(s) = F 1 (0)A(s) + F 2 (0)B(s), G(s) = G 1 (0)A(s) + G 2 (0)B(s). (62) 



The value of the form factors at zero is related to the derivatives of the pseu- 
doscalar masses with respect to the quark masses: 




- v 9m . 




(63) 

6',(0) = 1^1 =0 G.AO) (-^) • (64) 



m=0 



Up to now, we have followed the same line as Refs. [26,17]. But in these papers, 
the value of the scalar form factors at zero was derived using Standard %PT, 
i.e. supposing that the three-flavor quark condensate dominates the expansion 
of the pseudoscalar masses. We are going to allow a competition between 
the terms linear and quadratic in quark masses, so that the normalization 
of the form factors may become rather different from what is considered in 
Refs. [17,26]. In a similar way, the form factors that we will obtain could differ 
from the ones obtained by a matching with Standard one-loop expressions 
[27,30]. 

We consider here three models of T-matrix, proposed respectively by Oiler, 
Oset and Pelaez in Ref. [31], by Au, Morgan and Pennington in Ref. [32], 
and by Kaminski, Lesniak and Maillet in Ref. [33]. These models fit cor- 
rectly the available data in the scalar sector under 1.3 Gev, as discussed in 
Refs. [17,26,31]. However, they have to be corrected for very low and very high 
energies, as discussed in Ref. [17]: chiral symmetry constrains the nn phase 
shift near the threshold, and the asymptotic behavior of the phases shifts has 
to be changed to to insure existence and unicity for the solution of Eq. (61). 



5.3.2 Contribution of the first integral 

If we put Eq. (62) into Eq. (57), we obtain the spectral function as the sum 
of two contributions: 



Im II(s) =J W \ K 




32?r i=l,2 



<5 _ AM 2 

i -A i (s)B*(s)9(s-AMf) 



(65) 
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K 



Ml 
Ml 



16tt 



i=l,2 



AM} 



B i (s)B:(s)9(s-AM^ 



where the logarithmic derivatives of the masses are denoted: 



m s (dM\ 
Ml Vdm 



s / m=0 



m s dM 2 P 
Ml dm s ' 



1p 



m 
Ml 



'ml 

dm 



m=0 



(66) 



The two bracketed functions in Eq. (65) can be plotted: the first one is called 
"type AB*" , the second one "type BB*" . It is also interesting to study how 
these two contributions cancel each other inside the spectral function, by tak- 
ing the Standard tree-level estimates: 7^ = 1, Ax = 1 — M^/(2Mf C ) and 
= M%/(2Ml). A peak, corresponding to the narrow resonance /o(980) [8], 
arises with a height depending on the models: Ref. [31] leads to a smaller peak 
than Refs. [32] and [33]. 

The integral between and si in the sum rule Eq. (49) can be written, using 
Eq. (65): 

"1 - - . -2 



1 r If s \ M* 

- dslm U(s) - 1 = jAkIab + IkXk^Ibb, (67) 

■k J s V s / M± 

71 



where X X y — -M-xy ~ Mxy/s , involves the moments: 



M% = ^r 2 ]dss k Y. J S -^Ms)B*(s)9(s - AM?), (K 

*=1,2 



M * B= lkt I dS J S -^B^)BKs)6(s-AMr). (69) 



i=l-2 



Notice that we solve Omnes-Muskhelishvili equations to obtain the scalar form 
factors of the pion and the kaon in the limit m — > (and m s fixed at its 
physical value). But we consider T-matrix models fitting experimental data, 
with up and down quarks with their physical masses. The limit m — > sets the 
7T7r-threshold to zeroQ, changes tut phase shifts near the threshold and shifts 

5 For m — > 0, the cut along the real axis starts at s = 0. However, the integral 
Jq S ° ds (1 — s/so) ■ Im II(s)/s is convergent, since for s — » 0: Fi(s) — ► -Fi(O) and 
Gi(s) ~ Gi(0) • s, leading to: 

I m n(s)~ — — r ^^F 1 (0)G[ (0) -a. (70) 



25 



200 



100 







-100 







; 







.7 

) 
1 
i 









/ 
/ 

/ 1 
/ / 

/ / 
/ / 








1 


/ ^^^^^ / 
\ / 
\ / 








X y 





0.5 1 1.5 2 



50 



-50 











71 








\ N ^** 
\ N ^^"^ 

i \ 

\ ' 
i \ ' 
\ ' 
\ / 
\ / 


> ■ — 

/ 




\ / 
\ / 





0.5 1 1.5 2 



Fig. 6. Up: Contributions of type BB* (positive) and type AB* (negative) to 
the spectral function. Down: example of spectral function, obtained with 7^ = 1, 
X K = 1 - M%/(2M%) and 7^ = M%/(2M^). In both cases, we plot the results 
for the T-matrix models of Ref. [31] (solid lines), Ref. [32] (dotted) and Ref. [33] 
(dashed). 

slightly the KK threshold. Such modifications should not alter significantly 
the general shape of the spectral function. In particular, the integral of the 
spectral function, dominated by the / (980) peak, should be affected only 
marginally when T m ^ is considered instead of T m ^ . 
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5-4 Second sum rule : Si < s < sq 



The contribution of the integral below s\ is positive and dominated by the 
/o(980) peak. But according to Sec. 5.2, II is superconvergent, and the integral 
of the spectral function from to infinity vanishes. Im II(s) should therefore 
become negative in some range of energy. In particular, negative peaks should 
rather naturally appear in the spectral function, in relation with massive scalar 
resonances like / (1370) and / (1500) [8]. 

Let us suppose that the spectral function is negative for si < s < s .Q The 
contribution from the intermediate region in Eq. (48) can be estimated from: 

SO 

-J' <-- fds Im n(s) - (l - —) < —J', (71) 

S 7T J S V S Q J Si 

si 



where J' is the integral: 

SQ 

J' = - fds Im n(s) ( 1 - - ) , (72) 

7T J \ SqJ 

s\ 

which satisfies the sum rule: 

- Ids Im n(s) (l-—)+J' + — I ds n(s) (l - —) = (73) 
7r J V so J 2in J \ sr, J 

\s\=s 



The first integral in Eq. (73) can be computed from the spectral function 
obtained in the previous section: 

- J ds Im II(s) (l - -J = I^Ab + 1k\k^1' B b, (74) 



o 



with X' XY = M-xy ~ -Wxy/so involving the moments Eqs. (68)-(69). 

The contribution from the complex circle [third integral in Eq. (73)] can be 
estimated through Operator Product Expansion (OPE), using the method 
described in the following section: 



6 If the spectral function is partially positive in this range, our hypothesis will end 
up with an estimate for the second integral that will be smaller than its actual 
value. In that case, we would underestimate the difference X(2) — -X"(3). 
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— f ds U(s) ( 1 - — 
2m J w V s 

|s|=so 



(75) 



9[l-2C(3)]^X(2K( So )a 2 ( So ) 



Pa 7(7 + 1) /V 

2 " 8 3 



"2 U 



a 2 (s ) + 



9[l-2C(3)]^X(2K( So )a 2 ( So ) 
x [1 + 6.5 • a(s ) - 25.125 • a 2 (s )]- 



(76) 



5.5 High-energy contribution : \s\ — s 



We want the contribution of the integral on the large circle: 

/C = J_ / ds n(s) - (l - -) = — /d0 (1 + e* e ) n(p 2 = -s e*).(77) 
2m J s \ s / Ztt J 

\s\=S -7T 



The factor (1 — s/s ) suppresses the contribution stemming from the time-like 
region around So, so that we can use in this integral the Operator Product 
Expansion of II [34] . Once Renormalization Group Improvement is applied to 
Eq. (55), the QCD renormalization group invariant m{uu) gets the coefficient: 



a 2 (P 2 )mj(P 2 ) = a 2 (s )m 2 s (s ) x 



a(s ) 



8/60+2 



(78) 



with a(s) = a s (s)/n and b = 11 — 2Nf/3 = 9. The integral Eq. (77) becomes: 



9[1-2C(3)] Fl v ^m 2 (s ) 



2tt 



Ml 



X(2)^^-[a(s )]-^ f ds (l + e^ e ) a^(s e i9 ),(79) 

Sq J 



with 7 = 2 + 8/6o = 2 + 8/9. To compute this integral, we expand a(P 2 = 
s e' t9 ) in powers of a(s ). The behaviour of a(t) (t complex) depends on the (3 
function: 



-P[a(t)] = -p a 2 -fca 3 + 



7T 



33 - 2N f _ 9 
6 ~ 2' 



306 - 38^ 
A " 24 " 8 - 



(80) 
(81) 
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The expansion of a(s e i61 ) is: 



a(s e ie ) = a(s ) - ^Po9a 2 (s ) + 



a 3 (s ) + O(a 4 ). (82) 



We get: 



/C = 9[l-2C(3)]^X(2)^a 2 ( So ) 



X 



' Pol f X 



Pint ] 7(7 + 1) 



8 



= 9[l-2 C (3)]^X(2) ^ 



m^so) 2 



7T 



~2 0l 



(83) 



a 2 (s ) + • • • 



a 2 (s )[l - 6.5 • a(s ) + 48.236 • a 2 (s )]. 



This negative contribution is strongly suppressed by a 2 and m 2 /s . We have 
considered here m s ~ 200 MeV, but the contribution of this integral is so 
small that the error due to m s and a s can be neglected. Notice that duality 
is not supposed to arise in the scalar sector for as low energies as in other 
channels, due to a probably large contribution from the direct instantons in 
this sector [35]. 



6 Results 



6. 1 Logarithmic derivatives of pseudoscalar masses 



The logarithmic derivatives of the masses are obtained from the expansions of 
^ 2 M 2 : 



Fl and FlMl- 



m s (dM 2 \ m s dM 2 m (dM 2 \ 

AP Ml\dm s ) m=() Ml drn^ 1? Mj, { dm ) m=Q 1 } 



The corresponding expressions are given in App. C.l. We have A^ = since it 
is proportional to the derivative of M 2 with respect to m s in the limit m — > 0. 

We would obtain at the one-loop order in the Standard framework: 



7 '~ 1, 7 *~2§ = ' 04 ' 7 '~^ = - 02 - <85) 
Ml , 4MI , , 

A - ~ 1 " 2MJ = °' 96 - A " ~ M* = (86) 
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IK 
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0.876 


0.090 


0.081 
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0.078 


0.068 


0.958 


0.070 


0.059 


0.3 


0.920 


0.082 


0.071 


0.970 


0.069 


0.058 


0.995 
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0.049 


0.5 


0.946 


0.074 


0.062 


0.992 


0.060 


0.048 


1.015 


0.051 


0.038 


0.7 


0.967 


0.064 


0.052 


1.009 


0.050 


0.036 


1.029 


0.039 


0.025 


0.8 


0.975 


0.059 


0.046 


1.016 


0.044 


0.030 


1.034 


0.033 


0.019 


0.9 


- 


- 


- 


1.021 


0.038 


0.024 


1.038 


0.027 


0.013 


0. 


0.892 


0.080 


0.072 


0.943 


0.070 


0.060 


0.970 


0.063 


0.053 


0.3 


0.941 


0.072 


0.062 


0.990 


0.061 


0.050 


1.011 


0.053 


0.041 


0.5 


0.965 


0.063 


0.051 


1.009 


0.051 


0.038 


1.029 


0.042 


0.029 


0.7 


0.980 


0.052 


0.039 


1.020 


0.038 


0.025 


1.037 


0.029 


0.015 


0.8 


0.983 


0.046 


0.032 


1.022 


0.032 


0.018 


1.037 


0.022 


0.009 



Table 3 

Logarithmic derivatives with respect to m : 7^, 7^ and 7^, as functions of X(3) and 
r, for Fq=85 MeV (upper part) and Fq=75 MeV (lower part). 
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\k 
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A T; 
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1.438 


1.398 


1.455 


1.416 


1.467 


1.428 
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1.391 


1.346 


1.383 
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1.366 


1.321 


0.5 


1.326 


1.276 


1.284 


1.231 


1.226 


1.169 
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1.236 
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1.148 
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1.041 
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0.833 
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1.341 
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1.332 
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1.302 


1.263 


1.287 


1.247 


1.261 


1.220 
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1.224 


1.177 


1.163 


1.111 


1.086 


1.029 


0.7 


1.109 


1.053 


0.991 


0.931 


0.858 


0.796 


0.8 


1.041 


0.982 


0.894 


0.834 


0.739 


0.681 



Table 4 

Logarithmic derivatives with respect to m s : \k and X v , as functions of r and X(3), 
for Fq=85 MeV (upper part) and -Fq=75 MeV (lower part). 
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The logarithmic derivatives differ from these values because of the terms 
quadratic in the quark masses in the expansions of FpMp. The Tables 3-4 
collect values of these derivatives for r = 20, 25, 30, and F Q =75 MeV and 85 
MeV. We note that the values for F =85 MeV, r = 25, X(3) ~ 0.8 are in 
correct agreement with the Standard tree- level estimates Eqs. (85)- (86). 

We can notice that 7^ remains close to 1 if we change F , r and X{3). 7^ 
involves the linear m-dependence of the pion mass, which can be written as: 
F 2 M 2 = 2mS(2) + 0{m 2 ). Therefore, 7^ is related to the two-flavor quark 
condensate S(2), which is very weakly dependent on the values of F and 
X(3) [it is affected by the value of r, but our tables show only large values of r 
where X(2) does not strongly vary]. On the other hand, 7^ and 7^ are rather 
sensitive to X(3). These two logarithmic derivatives are 1 /r-suppressed when 
the three-flavor quark condensate is large. If X(3) decreases, these 2 quantities 
feel strongly the presence of large higher-order contribution. 

\ K and increase from 1 to 1.4 when X(3) decreases down to 0. If X(3) 
vanishes, the pseudoscalar masses are dominated by terms quadratic in m s . 
In that case, we would naively expect logarithmic derivatives to be twice as 
large as for X{3) ~ 1. To understand this discrepancy, it is useful to consider 
the second kind of logarithmic derivatives arising in Eq. (48): 

- m s dMp Mp n T ^ 

a * = m^7 = W Ap ' P = K ' n (8?) 



\ P and \ P are related, but the first is suppressed with respect to the latter 
by a factor Mp/Mp. For instance, in the limit of a vanishing three-flavor 
condensate, the expansion of FpMp begins with terms quadratic in m s , so 
that Ap is of order 2, whereas Ap = Ap • Mp/Mp is suppressed, and reaches 
lower values around 1.4-1.5. 

To obtain the values of Tables 3-4, we had to neglect the remainders of higher 
order ep and 5p. It is difficult to estimate the size of the resulting errors for 
the logarithmic derivatives 7p and Ap. Suppose for instance that ep/Fp and 
Sp/Mp are smaller than 10 %. To know the impact on 7p and Ap, we would 
have to calculate the values of the derivatives of ep and 5p with respect to 
m and m s . If we know only the size of e P /Fp and 5 P /Mp, it is hard to get 
any information about their derivatives, and to estimate the resulting errors 
on the logarithmic derivatives of the pseudoscalar masses. 
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Fig. 7. Sum rule : Range for -X"(3) as a function of r = m s /m for Fo=85 MeV, with 
the T-matrix models of Refs. [31] (up) and [32] (down). The results are plotted 
for ,si=1.2 GeV and so=l-5 GeV (solid lines), 1.6 GeV (dashed lines) and 1.7 GeV 
(dotted lines). The lines with white circles show the corresponding range for X{2). 



6.2 Estimate of X(3) and of LEC's 



Hence, two different estimates of X(2) — X(3) are available: the first one is 
the relation between X{2) and X(3) [Eq. (20)], the second one consists of the 
relation between X(2) - X(3) and 11(0) [Eq. (48)] and the sum rule for 11(0) 
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Fig. 8. Sum rule : Range for X(3) as a function of r = m s /m for i*o=85 MeV, 
with the T-matrix models of Ref. [33]. The results are plotted for si=1.2 GeV and 
so=1.5 GeV (solid lines), 1.6 GeV (dashed lines) and 1.7 GeV (dotted lines). The 
lines with white circles show the corresponding range for X{2). 

[Eq. (49)]. In both cases, the difference X(2) — X(3) can be expressed as a 
function of the observables and of F , r, X(3). This overdertermination can be 
viewed as a constraint fixing X(3) in terms of r and F , see Figs. 7-8. 

This analysis contains 3 sources of errors, i) First, we have neglected NNLO re- 
mainders in the expansions of pseudoscalar masses and decay constants. Their 
effect is easy to control in the relations between X{2) and X{3) [Eq. (20)] or 
F(2) and F(3) [Eq. (39)], but the situation gets more complicated for the sum 
rule Eq. (48) and for the logarithmic derivatives \p and jp. In the framework 
of Standard xPT, the authors of Ref. [6] noticed that the dependence on m s 
of S(2) is not really affected by two-loop effects. In addition, these effects 
have the same sign as one- loop contributions: if they were significant, they 
would increase (and not decrease) the gap between X{2) and X{3). A similar 
conclusion was drawn in Ref. [18]. The NNLO remainders are supposed here 
to be small, and they are not included in the results. 

%%) The evaluation of the sum rule Eq. (49) relies on an estimate of the integral 
between s\ and s . If we choose a couple (F ,r), we will not end up with one 
value for X(3), but rather a range of acceptable values that will also depend 
on the separators S\ < Sq. On the Figs. 7-8, the upper bound for X(3) remains 
stable for ^^sq > 1.5 GeV, whereas the lower bound depends strongly on s - 
When s increases, the lower bound of Eq. (71) is too loose to be saturated. 
A more stringent lower bound would be welcome. 
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r=ms/m 

Fig. 9. Sum rule : Low-energy constants ij=4...s(-^p) ■ 10 3 as functions of r = m s /m 
for Fo=85 MeV, s\=1.2 GeV and so=1.6 GeV, with the T-matrix models of 
Refs. [31] (up) and [32] (down). The values plotted on the left, along the verti- 
cal axis, are the Standard estimates stemming from Ref. [3]. 

in) The third source of error is the T-matrix used to build the spectral func- 
tion Eq. (65) for s < s±. Three different models of T-matrix have been used 
[32,33,31]. The central element is the shape of the /o(980) peak. Ref. [31] leads 
to the least pronounced effect. The two other models [32,33] lead to a higher 
/o(980) peak, a larger value for 11(0), and a smaller value for X(3). 
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Fig. 10. Sum rule : Low-energy constants Lj=4...8(M p ) • 10 3 as functions of r = m s /m 
for Fq=85 MeV, si=1.2 GeV and so=1.6 GeV, with the T-matrix models of 
Refs. [33]. The values plotted on the left, along the vertical axis, are the Standard 
estimates stemming from Ref. [3]. 

The range for X(3) is much narrower for large values of r, and can be even 
reduced to one value in the case of Ref. [31]. This range should be broadened 
if we took into account the errors related to higher orders in the expansion 
of pseudoscalar masses and decay constants. The value of F has no major 
influence on the constraint for [X(3),r]. For instance, choosing F =75 MeV 
would slightly shift the curves for X(3) towards the left of the graphs (r — > 
r — 2). Similarly, a change of y/s{ around 1.2 GeV does not affect strongly 
the results. If we choose v /sT=1.3 GeV, the convergence of the upper bound 
is slightly less good, but its values remain very close to Figs. 7-8. We should 
add a last comment for r ~ 25 (commonly used in the Standard framework): 
the values of X(3) correspond then to the half of X(2). We end up with a 
similar result to the one obtained in Refs. [17,18], but without relying on the 
hypothesis X(3) ~ I. 

The results of the sum rule for X(3) can be converted into bounds for Lj=4...8, 
plotted on Figs. 9-10 as functions of r, for = 1.2 GeV, — 1-6 GeV, and 
F = 85 MeV. For small r, the LEC's reach very large values: their definition 
from the low-energy behaviour of QCD correlators includes 1/B factors that 
make them diverge when X(3) — > 0. We notice also the large values of L 5 , 
L 7 and L 8 for r ~ 25, because the sum rule pushes L 6 (M p ) towards slightly 
positive values. The value of L 4 is not predicted by the sum rule: it depends 
essentially on the value fixed for F . For instance, choosing F =75 MeV would 
yield a slightly positive value for L 4 as r becomes large. 
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We have plotted on the left side of Figs. 9-10 the values of the LEC's of Ref. [3], 
which were derived assuming that X(3) is of order 1 and L 4 (M J? ) = Lq(M v ) = 
0. Let us remind that the values obtained for L 5 , L 7 , and L s are strongly 
dependent on these assumptions. The values of the LEC's of Ref. [3] hardly 
agree with the ones obtained from the sum rule, because the latter leads to 
positive values of L 6 (M P ) and to a small three- flavor condensate. 



6. 3 Slope of the strange scalar form factor of the pion 



Additional information about the decay constants is provided by the scalar 
form factors through a low-energy theorem. Consider the correlator: 

D%(p,q) = \\mjj e^-^(0\T{^(x) Aj(0) ss(yM {c K (88) 



where i,j = 1 . . . 3, and (c) denotes the connected part of (A 1 A{)(ss). The 
Ward identities p^D 1 ^ = r v D l ^ u = (with r = q — p) yield the Lorentz 
decomposition: 



D% jL = m skj{ K V 1 V i xVv ~ (p • f^VpJ'v + P 2 r ll r u - p 2 r 2 g^] 

+L[r^-(p-r)g^}}, (89) 

where K and L are scalar functions of p 2 , q 2 and r 2 . 

On the one hand, we have: 

d 



dm s 



i \ dV x (0\T{A;(x) Ai(0)}\0) 



(90) 



The correlator A^Al admits the following decomposition: 

i J d'xe^(0\T{A^x) Al(0)}\0) = 5^\p,p v - P 2 g^M P 2 ). (91) 

$(p 2 ) contains a pole at p 2 = stemming from one-pion states: 



= ¥- + ■■■ (92) 



where the dots denote contributions from the other states. We have therefore: 
D% (p, 0) = -8* [ P , Pu - g, vV 2 \ 1 1 + ...}. (93) 
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On the other hand, TJf u is dominated at low energy by the exchange of two 
pions between ss and each of the axial currents: 

D% = 25-F 2 (2) Gl ( q 2 )^- )+ ... (94) 



which contributes to K: 



K(p 2 , g 2 , r 2 ) = -2F 2 (2) f 2 l(g } + . . . (95) 



whereas L(p 2 ,q 2 ,r 2 ) receives no contribution. We compare Eqs. (93) and (95) 
for p, q — > to obtain: 

?m = 2FH2) Ita -^_^> = msG ; (0 ). (96) 



This low-energy theorem [17,18,26] provides a relation between the logarithmic 
derivative of F(2) with respect to m s , and the slope of the strange scalar form 
factor of the pion for a vanishing momentum. 

We can now exploit the solutions of Omnes-Muskhelishvili equations. Accord- 
ing to Eq. (62), we get for the slope of the form factor: 

DM 2 

m s G[(0) = m s G 2 (0)B[(0) = m s -^B[(0) = X K M 2 K B[(0). (97) 



B[(0) is computed by taking the derivative with respect to s at of Eq. (61): 



1 2 7 1 s'-4M 2 
^(0) = -E / dJ^wJ——LeV-AM?)B j ( i f). (98) 



The numerical resolution of Omnes-Muskhelishvili equations Eq. (61) yields 
the values of B(s) at the points of integration used for the Gauss-Legendre 
quadrature [17]. Hence, we can compute directly the integral Eq. (98) by the 
same integration method. 

On the other hand, Eq. (24) leads to: 

m s dF 2 {2) 1 [ ? rX (3)F 2 M%{ M 2 K M 2 K 



2F 2 {2) dm s F 2 



(99) 
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Fig. 11. Sum rule and slope of the strange form factor of the pion : Range for X(3) 
as a function of r = m s /m with the T-matrix models of Refs. [31] (up) and [32] 
(down). The results are plotted for s\=\.2 GeV and so=1.5 GeV (solid lines), 1.6 
GeV (dashed lines) and 1.7 GeV (dotted lines). The lines with white circles show 
the corresponding range for X(2). 

We see that Eq. (96) is an additional constraint, different from the sum rule 
Eq. (49). From the analysis of the pseudoscalar spectrum, we have concluded 
that all the quantities could be expressed (at the NLO) as functions of masses, 
decay constants, and 3 parameters F ,r, X(3). The sum rule was a first con- 
straint, fixing a range for X(3) depending on r and F = F(3). If we exploit 



38 




Fig. 12. Sum rule and slope of the strange form factor of the pion : Range for X(3) 
as a function of r = m s /m with the T-matrix models of Ref. [33]. The results are 
plotted for s\=1.2 GeV and so=1.5 GeV (solid lines), 1.6 GeV (dashed lines) and 
1.7 GeV (dotted lines). The lines with white circles show the corresponding range 
for X(2). 

the second constraint Eq. (96), we can obtain ranges for F(3) and X{3) as 
functions of r, plotted respectively in Figs. 11-12 and Figs. 13-14. These results 
can also be converted into values for the low-energy constants (see Figs. 15 
and 16). 

The values of X{3) are close to the ones obtained by the only application of the 
sum rule Eq. (49). The results obtained then for X(3) were not very sensitive 
to the valued chosen for F . We see also that the slope of the strange scalar 
form factor of the pion leads to rather small values for F(3) = F (around 70 
MeV) for r ~ 25. This result is in agreement with the small positive values 
obtained for L 4 (M p ). This increase of L 4 (with respect to the previous analysis) 
comes with a decrease of L 5 . In Ref. [26], the analysis of the same form factor 
led to d F = m s /F{2) ■ <91ogF(2) /dm s = 0.09. In the framework of Standard 
xPT, such a value corresponds to L 4 (M p ) ~ 0.4 • 10~ 3 , i.e. F(3) of order 75 
MeV. This question has also been discussed in Refs. [17,18,27]. 

However, these two constraints do not demand the same accuracy for the scalar 
form factors. The sum rule involves the integral of the spectral function Im II 
up to 1.2 GeV, which is dominated by the /o(980) peak. The global shape of 
the spectral function (and more precisely around 1 GeV) is the crucial element. 
For the low-energy theorem, we are interested in the slope of a form factor at 
zero, i.e. low-energy details. The resulting constraint may be less stable than 
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Fig. 13. Sum rule and slope of the strange form factor of the pion : Ranges for 
[F(3)/F 7T ] 2 (no symbol) and [F(2)/F 7T ] 2 (white circles) as functions of r = m s /m 
with the T-matrix models of Refs. [31] (up) and [32] (down). The results are plotted 
for ,si=1.2 GeV and a =1.5 GeV (solid lines), 1.6 GeV (dashed lines) and 1.7 GeV 
(dotted lines). 

the sum rule. It seemed therefore preferable to split the analysis in two parts: 
the first one dealing only with the sum rule, the second one exploiting both 
constraints at the same time. 
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Fig. 14. Sum rule and slope of the strange form factor of the pion : Ranges for 
[F(3)/F n ] 2 (no symbol) and [F(2)/F n ] 2 (white circles) as functions of r = m s /m 
with the T-matrix model of Refs. [33]. The results are plotted for s\=1.2 GeV and 
so=l-5 GeV (solid lines), 1.6 GeV (dashed lines) and 1.7 GeV (dotted lines). 

6.4 Scalar radius of the pion 



The scalar radius of the pion (r 2 )J can also be obtained from the scalar form 
factors of the pion, considered out of the chiral limit [i.e. with the physical 



masses m s and m 



(m u + m d )/2}\ 



F 1 (s) = F 1 (0) 



l+ l -(r 2 Y s s + c^s 2 + ... 




(100) 



If we project F on the two solutions A and B, we obtain: 



2\7T 



(r 2 ) 



^i(0) 



A> M + §fB[(0) 

lvl ir ITT 



(101) 



where a third kind of logarithmic derivatives is involved (considered out of the 
chiral limit): jp = (9[logMj,]/(9[logm]. App. C.2 collects their expressions in 
terms of the low-energy constants. 

We are interested in a quantity describing the non-strange pion form factor 
around the threshold. It should be possible to neglect the KK channel with 
no major change in the results. This point of view is supported by a numerical 
estimate: B[(0) / A[(0) ~ 0.1 and (M^/M 2 ) x (7^/77?) ~ 1/2. If we restricted 
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r=ms/m 

Fig. 15. Sum rule and slope of the strange form factor of the pion : Low-energy 
constants Li = ^_^{M p ) • 10 3 as functions of r = m s /m for ^0=85 MeV, si=1.2 GeV 
and so=l-6 GeV, with the T-matrix models of Refs. [31] (up) and [32] (down). 
The values plotted on the left, along the vertical axis, are the Standard estimates 
stemming from Ref. [3]. 

— # 

our analysis to the tttt channel, only the first term (the solution A) would 
appear on the right side of Eq. (101). The scalar radius of the pion would be 
independent of r, X(3) and F in that case. Actually, the second term on the 
right side of Eq. (101), related to the KK channel, is responsible for a weak 
dependence of (r 2 )1 on r, X(3), F . We can use the previous results, where 
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r=ms/m 

Fig. 16. Sum rule and slope of the strange form factor of the pion : Low-energy 
constants £j=4...8(-^p) ■ 10 3 as functions of r = m s /m for ^0=85 MeV, si=1.2 GeV 
and so=l-6 GeV, with the T-matrix model of Ref. [33]. The values plotted on the 
left, along the vertical axis, are the Standard estimates stemming from Ref. [3]. 

X(3) and F are functions of r, in order to study the range of variation for 
the pion scalar radius: 

0.537 - 0.588 fm 2 Oller-Oset-Pelaez Ref. [31], 

0.567 - 0.630 fm 2 Au-Morgan-Pennington Ref. [32], 
0.592 - 0.650 fm 2 Kaminski-Lesniak-Maillet Ref. [33], 

to be compared to the estimates: 0.6 ± 0.2 fm 2 [37], 0.55 ± 0.15 fm 2 [1], 0.55 
to 0.61 fm 2 [30], 0.57 to 0.61 fm 2 [38] and 0.61 ± 0.04 fm 2 [39]. Notice that 
the matching of Roy equations with Standard \PT [39] relies strongly on the 
value of the scalar radius of the pion. 

Information about the scalar radius of the pion could be seen as an additional 
constraint on our system, since (r 2 )J is related to 2L4 + L5. The situation 
is similar to dp'- this kind of constraint could rather easily be affected by 
higher-order corrections. We are also obliged to consider it out of the chiral 
limit m — > 0. It seems therefore wiser not to use this constraint, until a new 
analysis would treat less crudely NNLO remainders. 
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7 Conclusions 



The LEC's of the effective chiral Lagrangian should be determined as accu- 
rately as possible in order to know and understand the pattern of SB%S. These 
constants have generally been estimated from the expansion of Goldstone bo- 
son observables in powers of quark masses, supposing (1) a dominance of the 
quark condensate and (2) an agreement with the large- N c picture of QCD. But 
such determinations of the LEC's could be modified if quantum fluctuations 
turned out to be significant. A symptom of large quantum fluctuations could 
be seen in the large violation of the Zweig rule in the scalar channel and in 
large variations of chiral order parameters (e.g. the quark condensate) from 
N f = 2 to N f = 3. 

First we have studied how the relaxation of the Standard assumptions (1) and 
(2) could affect the determination of the LEC's. To reach this goal, we have 
studied the expansion in quark masses of the Goldstone boson masses and 
decay constants. We have truncated these expansions to keep the first two 
powers in quark masses and we have supposed that higher-order remainders 
[0(m 3 uark ) for FpMp and 0(m^ USiy . k ) for Fp] are small. These expansions can 
be written using "effective" scale-independent constants that combine chiral 
logarithms and LEC's. FpMp involves S(3), F 2 (3) and constants related to 
Lq, Li and L 8 , whereas Fp is expressed in terms of F 2 (3) and constants 
corresponding to L 4 and L 5 . 

We have not considered these expansions in one-loop Standard xPT, since the 
three-flavor quark condensate S(3) is not supposed to dominate the expansion 
of pseudoscalar masses. We have not worked either in tree-level Generalized 
xPT, since we have included chiral logarithms. These relations between LEC's 
and experimental quantities (masses and decay constants) can be inverted to 
express L i= 4... 8 [and therefore X(2) and F(2)] as functions of F(3), r = m s /m 
and X(3). We have then studied a possible competition between the first two 
orders of the quark mass expansions, by admitting large values for the ZR 
violating constants L 4 and L 6 (larger than the expected values on a basis of 
large- N c arguments). 

The variation of the Gell-Mann-Oakes-Renner ratio X from Nf = 2 to Nf = 3 
is governed by Lq. The equality X(2) = X(3) (saturation of the paramagnetic 
bound) is realized for L 6 (M P ) = -0.21 • 1(T 3 . A three-flavor GOR ratio X(3) 
much smaller than 1 could be obtained for two different reasons. On the one 
hand, the ratio of quark masses r may be smaller than 25 (r < 20), which 
leads to small values of X(2), and then of X(3) [due to the paramagnetic 
inequality X(2) > X(3)]. On the other hand, Lq(M p ) may be larger than the 
value —0.26 • 10~ 3 saturating the paramagnetic bound for X. A slight shift of 
L 6 (M p ) towards positive values leads to a significant decrease of X(3), whereas 
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A(2) remains almost constant and unsuppressed. A(2) could thus be of order 
1 and A(3) much smaller than 1, provided that r is large (r ~ 25) and the 
Zweig rule is strongly violated for the correlator defining L 6 . 

A similar analysis has been performed for the decay constants. L 4 tunes the 
difference between F 2 (2) and F 2 {3): the equality is obtained for L^{M P ) = 
-0.37 • 1(T 3 . If L 4 (M P ) is heading for positive values, F 2 (2) and F 2 (3) split 
swiftly. For r = 25, the saturation of both paramagnetic inequalities [for F 2 
and X] yields X{2) = A(3) = 0.9 and F{2) = F(3) = 87 MeV. This "ultra- 
Standard" scenario corresponds to the minimal values of L 4 and L 6 (no ZR 
violation). A slight drift towards positives values could lead to very different 
chiral structures of the vacuum for Nf = 2 and Nf = 3, corresponding to a 
significant role of quantum fluctuations in SB%S. 

The pseudoscalar spectrum (masses and decay constants) by itself does not 
contain enough information to pin down the size of these fluctuations. This 
effect can however be estimated from experimental data in the scalar channel, 
through a sum rule. The difference between A (2) and A (3) is related to the 
correlator II of two scalar densities uu and ss at vanishing momentum. 11(0) 
can be expressed in terms of a sum rule constituted of three distinct integrals, 
i) We compute the first one, involving the spectral function Im II up to energies 
around 1.2 GeV, by solving coupled Omnes-Muskhelishvili equations for the 
scalar form factors of the pion and the kaon. The solutions depend on the T- 
matrix model used to describe the interactions between tctt- and AA-channels, 
and on a normalization of the form factors related to the derivatives of M n 
and Mk with respect to m and m s . ii) The second integral corresponds to the 
contribution of the spectral function Im II between 1.2 and 1.6 GeV, where 
we cannot trust the two-channel approximation anymore. A second sum rule 
is used to estimate roughly this integral, iii) The third integral is performed 
on a large complex circle, with a large enough radius to rely on the Operator 
Product Expansion (OPE) of II. 

The most significant contribution stems from the first integral: the /o(980)- 
peak leads to a large value for 11(0), and therefore to an important splitting 
between A(2) and A(3). If we fix A(3), r and F(3), we know A(2) and the 
LEC's Lj =4 . 8 , using our previous analysis of the pseudoscalar spectrum. The 
derivatives of M w and Mk with respect to m and m s can then be directly 
computed, since they involve A(3), r, F(3) and LEC's. The sum rule Eq. (49) 
can therefore be seen as a constraint, giving A (3) as a function of r and 
F(3). Several sources of errors could affect this sum rule : the higher-order 
remainders in the expansions of FpMp and Fp, the rough estimate of the 
integral in the intermediate energy range, the T-matrix model. The three 
models considered here support nevertheless a large decrease of A (3) with 
respect to A (2), corresponding to positive values of L 6 (M P ). The size of the 
splitting between the quark condensates depends on the height of the / (980) 
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peak in the spectral function. In the particular case of "Standard" inputs 
r ~ 25, Fq = 85 MeV, the results of Ref. [17] are confirmed: X(3) can hardly 
reach more than one half of X(2) for the three considered models[]. 

The scalar form factors of the pion and the kaon can be exploited in several 
different ways. For instance, L 4 [i.e. F(3)} is related to the slope of the scalar 
form factor of the pion at zero. This second constraint may be used to fix 
X(3) and Fq as functions of r. If the conclusions for X(3) remain unchanged, 
positive values of L 4 (M P ) are obtained, leading to a significant decrease from 
F{2) to F(3) (20 to 30%). The Zweig rule would be violated strongly for L 4 
and Lq. However, this second constraint is sensitive to fine details of a form 
factor (slope at zero), whereas the sum rule depends on the general shape 
of the spectral function Im IT [and especially on the presence of a high peak 
corresponding to the /o(980) resonance]. The scalar radius of the pion has also 
been computed, in agreement with former estimates. 

A large decrease of the quark condensate from 2 to 3 flavors could be under- 
stood in terms of chiral phase transitions [13]. One of these transitions could 
be triggered by a vanishing quark condensate. If the corresponding critical 
value n CTit (N c ) turned out to be close to 2-3, we should expect significant vari- 
ations of the quark condensate with Nf in the vicinity of the critical point. 
Moreover, in terms of eigenmodes of the Dirac operator, the quark condensate 
can be interpreted as a density of eigenvalues, whereas L 6 corresponds to fluc- 
tuations of this density. Near the critical point where the first vanishes, the 
latter is expected to increase significantly. Let us remind that this scenario is 
only a possible explanation for a large difference between X(2) and X(3). The 
large value of ZR violating LEC's might be caused by another (and unrelated) 
mechanism. 

Forthcoming experiments [40] on tttt scattering should pin down the value of 
X(2), which is strongly correlated to r. If X{2) turned out to be close to 1, 
they could also measure low-energy constants of the SU(2) xSU(2) Lagrangian, 
Z3 and / 4 [1,41]. However, these experimental values could not be used to fix 
SU(3)xSU(3) LEC's without assumptions on the size of the ZR violating 
LEC's L 4 and L 6 [13]. It would be possible to constrain more tightly L 6 through 
a more sophisticated analysis of the sum rule including bounds on X{2) (or 
equivalently r). However, this remains a very indirect determination of the 
three-flavor condensate. Direct experimental tests are necessary to investigate 
closely the chiral structure of QCD vacuum for three massless quarks, and to 
understand the role of quantum fluctuations in the pattern of SB^S. 



We remind however that this result is barely consistent with the Standard hy- 
pothesis of a three- flavor condensate dominating the description of SBxS. 
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A Spectrum of pseudoscalar mesons 



A.l Decay constants 



The decay constants are [1,36] 



F 2 + 2m£ + 2(2m + m s )£ 
1 F 2 M 2 
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16tt 2 F 2 
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with the scale-independent low-energy constants: 
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and: 
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327T 2 
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(A.8) 



The higher-order contributions are denoted by S-zFp. The effective constants 
are related to F through the relations: 
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The decay constants fulfill the relation: 
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A. 2 Masses 



The pseudoscalar masses are [1,36]: 



F 2 M 2 = 2mS + {2mm s + 4m 2 )Z s + Am 2 A (A.12) 
F 4 M 4 

+^[X(3)] 2 L + FX, 

F 2 K M 2 K = (m, + m)£ + (m s + m)(m s + 2m)Z 5 + (m s + m) 2 A (A. 13) 



48 



1 FjM 4 

F4 



; f^(r + l)[X(3)] 2 L + F 2 K 5 K , 
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(A.14) 
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The higher-order remainders are denoted by 5Mp. The low-energy constants 
can be estimated using the relations: 
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with: 



A. 3 Pseudoscalar masses for m — > 

From the previous relations, one can derive low-energy constants from experi- 
mental data (pseudoscalar masses, F n and F K ) and 3 parameters: r, X(3) and 

F^.Fk F v ,m s £,m s i 
M n , M K , M v -> mfZ 5 , m^A, m*Z p 

In the chiral limit m — > 0, we will have to know the effective constants: 




(A.25) 



To compute Mp in this expression, we take the chiral limit of the mass expan- 
sions Eqs. (A.12)-(A.14). But these expansions involve the effective constants 
at m — 0, which leads to corrections containing logarithms of Mq/Mq: 

M 2 

M P = Y: a t X t + Y.D Q - log -2 (A.26) 

i Q M Q 

The equations Eq. (A.26) could be solved iteratively. Actually, Mq/Mq re- 
mains very close to 1. The calculation is simplified (and still accurate) if we 
compute in a slightly different way Mq in the logarithmic piece of Eq. (A.26). 
We start from Eq. (A.26), and we neglect the second (logarithmic) term: 

M Q = y £a i X i . (A.27) 



Mq is then directly computed from observables and F ,r, X(3). We put then 
these values of Mq in the logarithmic term of Eq. (A.26). We end up with 
values of Mp very close to the ones computed iteratively. These values will 
be used to compute the low-energy constants in the chiral limit Xj| m=0 using 
Eq. (A.25). 



Fjri Fxi F.q 

M K ,M V 



(A.24) 
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Fig. B.l. Self-energy diagram, leading to the same kind of integrals as in the OPE 
of II at the lowest order. 

B Operator Product Expansion for II 



Six integrals contribute to the Wilson coefficient of m s (uu) at the leading 
order in the strong coupling constant. The corresponding Feynman diagrams 
are drawn on Fig. 5. On each line, the left and right diagrams correspond to 
each other by crossing the gluonic lines. A simple change of variables in the 
integrals shows that the diagrams on the same line contribute identically to 
the Wilson coefficient. 

We want to consider the large-p 2 behavior of integrals like: 



J(M, {m },p) = I [g2 _ m ^[t-mi]^ (B '^ 



1 

x- 



[{k + q) 2 - ml] V3 [(p - q) 2 - ml]" 4 [(k +p) 2 - m|]^' 

These integrals are formally identical to the integrals arising in two-loop com- 
putations of self-energies, see Fig. B.l. 

The behavior of such integrals at large external momentum is known. The 
basic idea is to follow the flow of this large external momentum through the 
Feynman diagram, in order to Taylor expand correctly the propagators [24]. 
This procedure relies on the asymptotic expansion theorem [25] and can be 
formally expressed as: 

J r ~ X] Jr/-y o T {m . }i{(? . } J 7 . (B.2) 



7 



T denotes the whole graph, 7 are subgraphs into which the large external 
momentum may flow and T/j is the complementary graph of 7. For each 
subgraph 7 (see Fig. B.2), we write the corresponding Feynman integral J 7 . We 
perform then a Taylor expansion 7{ mi },{g 4 } with respect to the masses and the 
small momenta (external to 7 and not containing p). We combine the resulting 
"expanded" integral with the remaining graph r/7 and integrate over internal 
momenta. The asymptotic behavior of the whole integral Jr is obtained by 
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Fig. B.2. Subgraphs involved in the asymptotic expansion of the two- loop Feynman 
integrals. The solid lines constitute the subgraphs, the dashed lines correspond to 
the excluded propagators [24]. 

considering all the possible flows 7 for the large external momentum. 

We look for the leading order in 1/p 2 of the Wilson coefficient. All the subraphs 
do not contribute with the same power of 1/p 2 . In particular, the diagrams of 
type 5 do not appear in the Wilson coefficient of m s (uu) at the leading order 
in 1/p 2 . Gathering all the contributions, we obtain: 



n ^\r lS^i{[ 5 - 6 C(3)] + [5-6C(3)]-[l + 6C(3)]}, (B.3) 



where the bracketed terms correspond respectively to the contributions of the 
first, second and third lines in Fig. 5. The £(3) terms are related to subdia- 
grams of type 1 (see Fig. B.2), corresponding to two-loop massless integrals 



J(M,{o},p). 



C Logarithmic derivatives 



C. 1 Logarithmic derivatives for m — > 



To compute the logarithmic derivatives 7p and Ap: 
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we use the relations: 
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where X = \im m ^oX. 
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The logarithmic derivatives with respect to m are: 
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The logarithmic derivatives with respect to m s are: 

K = 0, (C.7) 
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C. 2 Logarithmic derivatives for m/0 



The same method can be used to compute the logarithmic derivatives involved 
in the scalar radius of the pion: 
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K 



Ml 



Ml 



4 4 
-4^ - -(4r + l)(r + 1)% - ^(16r 2 + lOr + 1)^ 



^ 2 M 2 X(3) 
F 2 F 2 64vr 2 



4 M 2 



3 "~° M 2 



d[F 2 5„ 



K 
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This linear system of three equations and three variables is easily solved to 
compute % and jk as functions of F , r and X(3). 
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